CERN-PH-TH/2006-185 
TPJU-10/2006 



A supersymmetric matrix model: 
III. Hidden SUSY in statistical systems 



G. Veneziano 

Theory Division, CERN, CH-1211 Geneva 23, Switzerland 

and 

College de France, 11 place M. Berthelot, 75005 Paris, France 

J. Wosiek 

M. Smoluchowski Institute of Physics, Jagellonian University 
Reymonta 4, 30-059 Cracow, Poland 

Abstract 

The Hamiltonian of a recently proposed supersymmetric matrix model has been shown 
to become block-diagonal in the large- A^, infinite 't Hooft coupling limit. We show 
that (most of) these blocks can be mapped into seemingly non-supersymmetric (1 + 1)- 
dimensional statistical systems, thus implying non-trivial (and apparently yet-unknown) 
relations within their spectra. Furthermore, the ground states of XXZ-chains with an 
odd number of sites and asymmetry parameter A = —1/2, objects of the much-discussed 
Razumov-Stroganov conjectures, turn out to be just the strong-coupling supersymmetric 
vacua of our matrix model. 
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1 Introduction 



In a recent series of papers we have introduced a supersymmetric quantum mechan- 

ical matrix model and studied its rather intriguing properties in the planar approximation, 
or, formally, in the large- A^, fixed-A limit (here is the size of our bosonic and fermionic 
matrices and A = g^N is the usual 't Hooft coupling) [5] . In particular, the model ex- 
hibits, at a critical value of A, Ac = 1, a discontinuous phase transition characterized 
by the emergence of new supersymmetric vacua on the strong-coupling side of the phase 
transition, and a consequent jump of Witten's index |6] across Ac- 

This property was first observed [H [2] in the lowest fermion- number sector of the 
model, F = (-F being an exactly conserved quantum number for all values of A) where 
one new supersymmetric vacuum emerges at A > 1. It was later realized that a similar 
phenomenon also occurs at F = 2 ^ , with two new supersymmetric vacua popping up at 
A > 1. In that same paper a deeper understanding of this unexpected feature was gained 
by considering the A ^ oo limit of the model. In this limit the Hamiltonian becomes 
block-diagonal in both F and boson number B, with blocks of finite size Af{F, B). 

Furthermore, by combining the strong-coupling limit with some combinatorics argu- 
ments [3], it was conjectured that the pattern found at F = 2 should generalize to all 
even values of F\ two new supersymmetric vacua would occur in each one of these sectors 
at large A. By going to infinite A, and by computing suitable supertraces [3], we can also 
guess which blocks should "hide" (for a given even F) the two new ground states: those 
with B = F ± 1, a conjecture confirmed by many numerical checks. Finally, once the 
infinite-A ground states are identified, their expression at finite A can be reconstructed 
through an explicit formal operation [1] , which is expected to lead to a normalizable state 
if and only if A > Ac. 

In this paper we consider again the A — oo limit, albeit for a different purpose. We 
will show that some of the finite blocks (including those where the new supersymmetric 
vacua occur) can be mapped into (seemingly non-supersymmetric) one-dimensional sta- 
tistical mechanics models with a finite number of sites and a periodic (cyclic) structure. 
More explicitly, we will map some sectors of our model into the XXZ Heisenberg chain 
with asymmetry parameter A = ±1/2. Moreover, it will be argued, and demonstrated 
numerically, that our system is also equivalent to a lattice gas of g-bosons in the limit 
where the quantum-deformation parameter q goes to infinity. 

Such kinds of equivalences are not new, the Thirring-sine- Gordon connection [7j being 
a famous example; however, the equivalence discussed in this paper has two novel features: 

• It connects a supersymmetric system with a (seemingly) non-supersymmetric one, 
hopefully revealing hidden supersymmetric features of the latter model. As an 
example, the ground state of the XXZ model with A = —1/2 is known to have 
amazing (partly proved, partly conjectured) properties [B] that could possibly be 
explained after realizing that such a ground state is just a supersymmetric vacuum. 

• It relates a rather abstract quantum mechanical matrix model, in the planar ap- 
proximation, to some well known statistical system in one-space one-time, therefore 
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providing an a priori unexpected physical (string-like?) interpretation of the former. 

In the next section we recall the definition of our model, its main physical properties, 
and its large- A limit. The equivalence with the XXZ chain and with the g-bosonic gas is 
discussed in the following two sections. We will end with a summary and a discussion of 
the possible consequences of this equivalence for the latter two systems. 

2 A planar supersymmetric matrix model 

The model is simply the iV — > cxd (planar) limit of an iV x iV matrix system defined by 
the following supersymmetric charges and Hamiltonian: 

Q = TT[fa^{l + ga% = Tilf^ {1 + ga)a], (1) 

H = {Q\Q} = Hb + Hf, (2) 
Hb = Tr[a^a + g{a^^a + a^a^) + g'^a^'^a'^] , (3) 

Hp = Tr[/t/ + g{f^f{a^ + a) + f\a^ + a)f) 

+ g\f^afa^ + f^aa^f + f^fa^a + f^a^fa)] , (4) 

where bosonic and fermionic destruction and creation operators satisfy 

[aij,a\i\ = 5ii5jk \ {fijfli} = ^u^jk] i,j,k,l = l,...N, (5) 

all other (anti) commutators being zero. 

Models of this type result from the dimensional reduction of D = (1 + l)-dimensional 
gauge theories. For example, two-dimensional Yang-Mills gluodynamics, when reduced 
to QM, is described by a free Hamiltonian Hyuz = — Tr[(a — a^)^]/2 acting on the gauge- 
invariant subspace of the whole Hilbert space [S] . The Hamiltonian ([2H1]) was designed to 
illustrate a new, general method [T] of finding the spectrum of gauge systems at infinite 
number of colours A^. It turned out, however, to have an interest of its own, by exhibiting 
the following non-trivial properties: 

• Since (jl]) conserves the fermionic number F = Tr [/'''/], the system can be studied 
separately for each F. 

• The planar model is exactly soluble in the F = 0,1 sectors, i.e. the complete energy 
spectrum and the eigenstates are available in analytic form. 

• There is a discontinuous phase transition in the 't Hooft coupling at A = Ac = 1. 
At this point the otherwise discrete spectrum loses its energy gap and becomes 
continuous. 

• An exact duality between the strong and weak couphng phases occurs in the F = 0, 1 
sectors. 
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• The system exhibits unbroken supersymmetry, i.e. there are exact, SUSY-induced 
degeneracies between bosonic (even F) and fermionic (odd F) eigenenergies. 

• In the weak couphng phase, A < 1, there exists only one (unpaired) SUSY vacuum. 
It lies in the F = sector and it is nothing else than the empty Fock state |0). For 
A > 1, however, the structure is much less trivial and more interesting: there are 
two SUSY vacua in each bosonic sector of the model. For F = the empty Fock 
state continues to be a null eigenstate, but it is joined by another, non-trivial, ana- 
lytically known ground state. For higher even F, two new non-trivial vacua appear. 
This is possible thanks to the rather intriguing rearrangement of the members of 
supermultiplets that occurs across the phase-transition point. 

We have established all these points for the F = 0, 1, 2, 3 sectors and believe that this 
structure persists for arbitrary F. This expectation is borne out by considering the infinite 
A — i> oo limit of the Hamiltonian (jlj) . Since this is also the limit in which our Hamiltonian 
reveals the above-mentioned connections to statistical mechanics, let us recall the strong 
coupling limit of our system [1] in a little more detail. 

Define the (appropriately rescaled) strong coupling SUSY charges by: 

Qsc = lim ^ Q = ^Tr(/at') , = ^Trifa') . (6) 



The corresponding strong-coupling Hamiltonian is just their anticommutator. Doing the 
algebra carefully, and throwing away terms that do not contribute in the large- limit, 
we find: 

Hsc = lim = ^Ma^'a^ + fafa^ + faa^f + ffa^a + fa^fa] . (7) 

This can be further simplified with the aid of the "planar calculus" rules derived in [1], H] ; 
namely, the third and the fourth terms must be brought to normal form, giving: 

fijajkahifh = flj{a\iajk + 5ji6kk)fii iVTr[/V], (8) 
fijfjkahiaii = fljO^kifjkaii . (9) 

In the above relations we have neglected terms in which the normal ordering does not 
match the trace-ordering, since such structures do not contribute to leading order in 
[HUj. In conclusion, the strong-coupling Hamiltonian reads0 

i^sc = Tr[/t/ + l(at'a2 + a^faf + fa^fa)]. (10) 

Remarkably, (fTOl) conserves both F and B = Tr[a^"a]. As such, the infinite Hamiltonian 
matrix splits into finite blocks labelled by the number of fermionic and bosonic quanta. 



^To be precise H^c contains also a Tr(/t)Tr(/) term. It is subleading except when it acts on 
\B = 0, F = 1) and plays an important role for assuring degeneracy with |_B = 2, F = 0). 
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{F, B). In each such sector the leading-order (i.e. planar) basis is generated by the single 
trace of a product of elementary creation operators. We may thus represent the generic 
state in a block of given B and F in the form: 

K,n,) = -^TT[a^'^\fra^'^\fr\..{f^r'^]\0)-m,,n. > 0. (11) 

This is a state with B = ^mi bosons and F = ^ni fermions, Mn being a normalization 
factor. Owing to the Pauli principle not all configurations of {mj,nj} define a legitimate 
state. The detailed rules for counting such states (Pauli-allowed necklaces) follow from 
Polya's theory and have been developed in [3] : they give the dimensionality of each block of 
the strong-coupling Hamiltonian ffTOl) . Table 1 shows a map of the first few (F, B) sectors 
together with their sizes. Since, in the strong-coupling limit, the supersymmetric charges 
connect states with the same value oi B + 2F, taking supertraces at fixed B + 2F is a way 
to check whether all states with such a value of i? + 2F are paired into supermultiplets 
or not. Such an analysis was carried out in [3] and revealed a bosonic excess by one unit 
for (and only for) B + 2F = ±l(mod 6). 

Using the planar rules developed in [H |2l H] the finite-dimensional Hamiltonian matrix 
can be readily calculated in each sector. Proceeding in this way, we have discovered 
"experimentally" the existence of a "magic staircase" - a distinct set of sectors (labelled 
by a bold face in Table 1) - where the zero-energy eigenstates are located!^. Indeed the 
magic sectors appear to lie at 

B = F±l,F = 2n^B + 2F = 3F±l = 6n±l, (12) 

in agreement with the supertrace considerations of [3] . This observation also explains the 
structure of SUSY vacua in the whole strong-coupling phase. Indeed, we were able [3] to 
give a (formal) way to express SUSY vacua at finite A in terms of those at A = oo. The 
connection is expected to lead to normalizable states only at A > 1. 

It remains, however, to understand why the strong SUSY vacua exist solely in the 
magic sectors f|T2l) . This puzzle finds its solution upon mapping our system into the XXZ 
Heisenberg chain. 



3 Equivalence with the XXZ Heisenberg chain 
3.1 The XXZ model 

The XXZ system [TOl [TTl [T2] is a one-dimensional periodic lattice of size L with a spin 
1/2 variable residing on its sites. Its Hamiltonian depends on an "asymmetry" parameter 

^The dimensionality of these sectors turns out to be given by Catalan's numbers [3]. 
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Table 1: Sizes of gauge-invariant bases in the {F,B) sectors 
A and can be written in terms of Pauli's matrices (with o"^ = |(cr^' ± ia"^)) as: 

i=l 

= - E + ^r<i + Y^f^i) = - (o± + f o.) , (13) 

where the site L + 1 is identified with the site 1. A convenient basis of states, labelled by 
two sets of integers {mi,ni), is the following: 

|m„n,) = |(0)"^Hir---(0)'"'-(l)"'-) , r>l,m„n,>0, EK + ^«)=^' (14) 

where a (0)™ (resp. a (1)"") indicates a sequence of m (resp. n) spins pointing down 
(resp. up). The state is defined up to cyclic permutations and the permutation giving 
the smallest binary number can be taken as representative. To the set fll4p one still has 
to add two states: (0)-^ and (1)^. 

The action of Oz on such states is very simple: 

Oz \mi, Hi) = (rieq - riopp) \mi, rii) , (15) 

where rieq (^opp) is the number of equal(opposite)-spin nearest neighbours. It is easy to 
see that rieq = L — 2r while riopp = 2r. Therefore: 

Oz\mi,ni) = {L - Ar) |mj,nj) . (16) 
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The action of 0± is only a bit more complicated: 

0±\mi,n^) = ^\m'-,n'-) , (17) 

where the sum extends over the sets {m'j,nj) that are obtained from {mi,ni) by inter- 
changing, in turn one by one, any pair of opposite-spin neighbours. As a result, the 
Hamiltonian commutes with the z-component of the total spin and is block-diagonal with 
blocks of given m = ^ mj and n = ^ rij (n + m = L). In conclusion: 

^xxz l^^i^ni) = -^(L-4r) \mi,ni) - ^\m'j,n'j) 

= L - 2A (L - r)^ \mi, m) - ^ \m'j,n'^) . (18) 

It is useful to define a rescaled XXZ Hamiltonian by: 



4xz = -^4xz - + 



so that: 



i=i ^ ^ 

For the special values of the asymmetry parameter A = ±i we get: 

m^' = -lL + Hi^K (21) 

wher^: 

i/W \mi, rii) = {L-r) \mi, n^) ± ^ |m;, n;.) . (22) 

We will now argue that the strong-coupling Hamiltonian of our supersymmetric matrix 
model can be mapped into if*^^) for some specific values of 5, F that include those of the 
magic staircase. 

3.2 Proof of the equivalence 

Let us start by splitting the strong-coupling Hamiltonian fITOl) as follows: 

H,c = + i/g , (23) 



We are grateful to Don Zagier for having pointed out to us this simple version of the XXZ Hamiltonian 
and for having suggested a possible relation to our matrix model. 
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where: 

HSi = Tr(/V) + ^Tr(at^a2) , (24) 

^sc = ^ + a^Paf] , (25) 

and the labels d (o) stand for diagonal (off-diagonal) pieces of the Hamiltonian. 

Let us now consider the action of Hsc on the generic state flTT]) whose similarity with 
the XXZ states ( IT^ is evident. The planar rules for doing that were discussed in [1], [2]. 
For H^'^} we simply find [F = n): 

Hf^\ni,mi) = ^F + ^(mi-l)j |ni,mi) = (L-r) |ni,mi). (26) 

The action of -f^gc? on the other hand, is precisely to interchange each fermion-boson pair 
of neighbours, i.e. an action very close to that of 0± in the XXZ model. However, since 
sign problems can arise, we have to treat various cases separately: 

• Fodd 

In this case there are no relative signs from different cyclic orderings in (fTTj) and 
therefore the action of ifg^ on them is exactly the same as that of 0± on the states 
f fT^ . The relative sign of ifg^ and ifg^ in ifgc concides with the one in H^+\ We 
thus obtain: 

H,o ^ H^^^ = + \l = -i^&f ^ + \l . (27) 

• F even and i? odd 

In this case different cyclic orderings in ffTTj) do carry relative signs. We will argue 
that, by a suitable choice of the basis, we can bring ifgc to agree, up to a shift, with 

o-(-l/2). 

-"xxz ■ 

^ H(-^ = ) + h . (28) 

The argument goes as follows: for the XXZ model let us choose as basis of states: 

(nil . . . 1100 . . . 000), (nil . . . loio . . . ooo), . . . (29) 

The corresponding basis for Hsc is then taken to be: 

l/VW^ . . . /V^«^at . . . atatat), -|/t/t/t/t . . . /tat/tat . . . a^a^a^), . . . , (30) 

i.e. we put in correspondence the ones (zeroes) in (!29l) with the fermions (bosons) 
in ( l30ll . while assigning a sign (— I)'' to a state in the latter set if it is obtained from 
the first state by interchanging k boson-fermion pairs. 

It is quite obvious that, on this convenient basis, H^c gives the same matrix elements 
as H^~\ modulo the possibility that ifgc produces a cyclic permutation of a state in 
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the above basis. However, even in this case, the correspondence works fine thanks to 
the fact that the relative sign (—1)^ originating from Fermi statistics (where p is the 
number of fermions to be interchanged) is equal to the relative sign (—1)'^ counting 
the number of boson-fermion interchanges. This is so since taking a fermion from 
the last to the first entry corresponds to an odd number of fermion-fermion and 
boson-fermion interchanges, while doing the same with a boson involves an even 
number of interchanges of each type. Note that for this to be true it is essential 
that B be odd and F be even, and, therefore, that L be odd. The importance of L 
being odd was much emphasized in [8]. 

• F and B even 

In this case we have found no simple relation between H^c and Hxxz for any choice 
of A. 

• We may add here a side remark: in the case of F odd, if we define a new (non- 
super symmetric) Hamiltonian: 

Hsc = H^sc - 4i , (31) 

we also find: 

Hsc ^ H("^ = ) + . (32) 



3.3 Implications of supersymmetry on the XXZ model 

The final form of our equivalences reads: 



-4xf ^ + lL, F odd 

r(-l/2) , 3 
'XXZ ' 4 

Eq. (IT^ implies that the parameters of both systems are related as follows: 



L = F + B, (34) 
S. = l{F-B), (35) 

where 5*^ is the conserved component of the total spin. In addition, the spectrum on the 
spin side should be computed in the sector that is invariant under the lattice shifts. 

We checked Eqs. fl33l) for all sectors with 5 < F+B < 9. Everything works as expected, 
including the magic sectors with a zero eigenvalue. Since the spectrum of Hsc is positive 
semi-definite, the 2nd of Eqs. (l33l) gives a simple and elegant proof that, for L odd, the 
states considered in are indeed ground states of the XXZ chain. The fact that they 
have Sz = ±| just corresponds to our SUSY ground states having B = F±l and F even. 
Moreover, for sectors with F and B even, where the equivalence is not expected to work, 
we indeed find different spectra in the two models. 
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Supersymmetry also implies that the (non-vanishing) spectrum of Hgc in the sector 
{F, B) has to be contained in the spectra of the "neighbouring" sectors with (F=i=l, S±2). 
If we take B odd, i? ± 2 is also odd and therefore we are always in a situation in which 
we are able to connect Hsc{F,B) to Hxxz{n,m). Obviously, n = F and m = B are 
conserved by i^xxz as well as by Hsc- After some trivial algebra, we get the following 
predictions (with integer /i and v ensuring even/odd m and n respectively): 

1. The spectrum of H^^l^\2v+1, 2/i+l) is contained in the combined spectra 

of -H^~^l^\2u, 2/i + 3) - 3/4 and -H^^^l^\2u + 2, 2^ - 1) + 3/4, 
and: 

2. The excited spectrum of i7^x2^^(2z/, 2/i + 1) is contained in the combined 
spectra of -H^^^l^\2u - 1, 2/i + 3) + 3/4 and -H^^^l^\2u + 1, 2/i - 1) - 3/4. 

Since the Hxxz{n,m) Hamiltonian is symmetric under the exchange of m and n, all 
variations of the above relations resulting from the interchange m ^ n are also valid. 
We made extensive numerical checks of these (to our knowledge novel) relations between 
different XXZ models at different values of A. 

3.4 Bethe ansatz solutions for the lower sectors 
of the XXZ chain 

The XXZ chain is integrable [12]. In particular the eigenenergies of ifxxz(A) are given 
exactly by the Bethe ansatz [TOl [TBJ : 

A 

^xxz(A) = -L- + 2mA - 2 J] cosp,, (36) 
where the momenta — vr < pj < n satisfy the following set of Bethe equations: 

With m denoting the number of down spins in a chain. For the supersymmetric model 
this translates into 

Esc = F + 2 J2L cospj, for F odd , A = +^ (38) 

2 

Esc = F — 2 '^f^i cospj, for F even, and B odd , A = — . (39) 

2 

Given the Bethe momenta pi, all eigenvectors can also be explicitly constructed. The 
whole problem reduces therefore to the solution of the non-linear equations (1371) . Con- 
sequently, the existence of the magic staircase with its supersymmetric vacua follows 
directly from the Bethe ansatz solution of the XXZ chain. The literature on the latter 
subject is huge, see e.g. [I2]-[T6] for more references. We shall content ourselves here 
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with formulating only an "ansatz within the Bethe ansatz" , which reduces the number of 
Bethe momenta needed to find our SUSY vacua. Solving numerically Eqs. (1S7I ) for a few 
low-F magic sectors, we have found that the zero-energy momenta satisfy (see Table 2): 

Pi = 0, p2k+i = -P2k, k = l,...,{B-l)/2. (40) 

In words, there is always one zero momentum, and the remaining ones come in pairs 
with opposite sigij^. We conjecture that this configuration gives the zero-energy state for 
arbitrary B = F ±1 and even F. 





Af{F,B) 




P2/vr 


Pi/n pe/n 


(4,3) 


5 


0.0 


1/3 




(4,5) 


14 


0.0 


0.260669 


0.558585 


(6,5) 


42 


0.0 


0.210767 


0.432222 


(6,7) 


132 


0.0 


0.178899 


0.364275 0.583645 



Table 2: Bethe momenta of the supersymmetric vacua in some of the lowest magic [F, B) 
sectors, see Table 1 

It is perhaps amusing that the ansatz fj40l) allows the Bethe phase factors to be derived 
for the first three vacua listed in Table 2 in analytic form. Assuming that they correspond 
to zero-energy eigenstates, we look for simultaneous solutions of Eqs. ( 137|) and E{F, B) = 
0. This problem can be solved algebraically. Defining Xi = e*^% we find: 

F = 4, B = 3 : 
F = 4, B = 5 : 




"^Recall that the magic sectors only occur for odd values of B. 
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corresponding to the algebraic representations of the phase factors given in Table 2. Al- 
though it is not self-evident, it can be proved algebraically that the above numbers have 
modulus 1. 



3.5 The Razumov— Stroganov states and supersymmetry 

Interestingly, the magic staircase, with its supersymmetric vacua, connects directly to 
properties of one-dimensional spin chains. Beginning with the classic paper of Baxter 
[llj, some simple eigenvalues of the XXZ Heisenberg chain were discovered. In particular. 



to any finite, odd L and made several conjectures on the properties of the eigenvector 
that corresponds to the above-mentioned ground state (for recent developments see, e.g. 



It follows directly from fl33l) and fl35|) that these states are nothing but the super- 
symmetric vacua of our planar model with B = F ± 1,F even. Hence Razumov and 
Stroganov's above-mentioned result guarantees the existence of one bosonic SUSY vac- 
uum in each one of our magic sectors. Hopefully, this hidden supersymmetry will help 
understanding and/or proving the other amazing (and so far mostly conjectured) proper- 
ties of Baxter's ground states. 

It is tempting to say that the two families of vacua, i.e. those with some given even F 
and B = F ± 1, are related by the usual inversion, cTj —Ui-, symmetry. However, this 
transformation has to be coupled with a change in the lattice size L — L + 2. This brings 
the issue of whether a change of L should be considered as a symmetry. At first sight this 
looks a little premature, and in fact, for the above vacuum solutions, it is pure semantics. 
However, this is no longer the case when we consider the implications of supersymmetry 
on the whole spectrum, including all excited states. The supersymmery generators act 
as F — »• F ± 1 and B —>■ 5 =p 2, corresponding to L — > L =p 1. Therefore they do relate 
the spectra of excited states on different lattices: the XXZ chain turns out to have a 
hidden supersymmetry with different members of its dynamical supermultiplets living on 
different lattices! 

4 Equivalence with a g-bosonic gas 

Surprisingly, our infinite-coupling planar system is exactly equivalent to yet another, 
and apparently entirely different, model. To expose this equivalence we use the original 
labelling of planar states with F fermions 




I15J). 




(41) 



Hsc = Tr[/t/ + ha^\' + a^faf + fa^fa)] , 



(42) 
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on a state (HTl) . Using the planar rules developed in it is easy to show that the first 

two terms do not change the initial state and give rise to the following diagonal elements: 



F 



(ni,n2, . . . ,nF\Hsc\ni,n2, . . . .np) = F + '^{n^ - 1 + 5„^,o) = B + ^5„,,o • (43) 

i=l i=l 

Using the same planar rules, the remaining matrix elements can be easily obtained from 
the action of the last two terms, e.g. 

Tr[aV^a/]|?^i,^2, . . . .np) = - 1, ris, . . . , rij^ + 1) + ^^I'^i + 1, - 1, • • • , '^f) + 

-jj^\ni,n2 + l,n3 - 1, . . . + . . . + -^|ni,n2, . . ..np-i + l,nF - 1) , 

(44) 

i.e. this operator annihilates one bosonic quantum in a group i and adds one at i — 1, 
meaning at the cyclic left of = 1, Similarly the last term moves one quantum 

to the cyclic right of i. Here, Mn {-f^nf) are the normalization factors of the initial (final) 
states of our basis. They contain some powers of and degeneracy factors ds- When 
calculating matrix elements of the Hamiltonian Hsc, all N factors cancel with the in 
( H2|) and we are left only with the square roots of the ratios of corresponding dg factors. 

It turns out that this Hamiltonian also describes the following system. Consider a 
one-dimensional, periodic lattice of length F. Put at each lattice site a bosonic degree of 
freedom described by the creation/annihilation (c/a) operators ai, i = 1, F and use the 
harmonic oscillator basis. The states |ni, n2, rip) are described by the configuration of 
F integer occupation numbers as before. 

The new Hamiltonian reads 

F F 
j=l i=l 

where Ni is the usual operator of the number of quanta at a site i, Ni = a^itti. In words: 
the second term counts the number Uzer of empty (unoccupied) sites in a given basis 
state and returns this state multiplied by Uzer- B is the total number of bosonic quanta, 
B = ni + n2 + ... +np. The h\ (6j) operators create (annihilate) one quantum without the 
usual ^/n factors. Omitting momentarily the index i: 

6t|^) = |^ + l)^ h\n) = \n-l), 6|0)=0. (46) 

In terms of the usual a, operators they read: 

1 1 
= , b = g, and 610) = , (47) 

VnTT VnTT 

where again = a^a. The b operators have non-standard commutation relations: 

[6,61=5^,0- (48) 



13 



This Hamiltonian conserves the total bosonic number, as before. It is also invariant 
under lattice shifts and, consequently, commutes with the lattice-shift operator U defined 
as: 



U\ni,n2, ...,nF) = \n2,n3, ...,np,ni .) 



(49) 



Therefore, the Hilbert space of states with fixed B can be further split into sectors with 
fixed eigenvalues of U: 



We claim that the spectrum of the above H, in the sector with Xjj = (—1)^"^, exactly 
coincides with our spectrum of Hsc, for all F and B. 

The main steps in understanding this equivalence are as follows: 

• Our planar states (HTl) are defined modulo Zp shifts. Without fermionic degrees of 
freedom this would be taken care of by requiring Xu = 1 for the bosonic system. 
The minus sign is the consequence of the fermionic operators in (HTj) : under the Zp 
shifts they acquire the phase (—1)'^"^. 

• For even F, the projection for the Xu = —1 sector plays another important role. 
Namely, it removes states that are not allowed by the Pauli principle. 

• Finally, the degeneracy factors required in fl44p are correctly taken into account by 
the linear combinations corresponding to the condition Xu = (—1)'^"^. 

The last point is particularly non-trivial. We have therefore double-checked this equiv- 
alence by diagonalizing both Hamiltonians in a range of sectors :3<F<7, 3<-B<7. 
All spectra are indeed identical, including again the supersymmetric vacua in the magic 
sectors. 

Notice that this second equivalence also works for the "bad" sectors with both F and 
B even. These are the only sectors where the Pauli principle is effective and eliminates 
some of the planar states. This is why there is no XXZ equivalence in these sectors. 
Nevertheless, because of the Xu = —1 projection, the bosonic- model equivalence applies 
to these cases as well. 

Interestingly, the system of "funny" b and operators turns out to be a particular limit 
of the so-called g-deformed harmonic oscillator algebra, well known in the literature [T71 
[TS] . The transitions (14^ (without the square roots) are referred to as asszsted transitions. 

The g-boson operators satisfy (for one degree of freedom) 



with = a^a, the usual occupation number operator. The b, 6^ operators are related to 
the standard a,a'^ by 




= 6^"^^, m = l,2 



(50) 



bb^ - q~%^b = 1, [N, b] = -b, [N, b^] = b^ , 



(51) 




(52) 
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where: 



l-q 



-2x 




—2 



(53) 



An alternative form of the commutation relations reads: 



[b,b^]=q 



(54) 



It can readily be checked that the usual harmonic oscillator algebra is recovered in the 
limit g — i> 1, 6, 6^ — s> a, a^. On the other hand, in the limit q ^ oo, b and b^ become our 
c/a operators satisfying (H6ll - (H8l) . 

An important point is that the g-bosonic Hamiltonian without the commutator (or 5) 
term in fH5|) is exactly soluble for all values of the deformation parameter q. However, 
with the additional 6 term, it is not. On the other hand, given the present chain of 
equivalences, we see that the above non-linear system oi q = oo-bosons is soluble in terms 
of the Bethe ansatz for the XXZ chain. Finally, and similarly to the latter case, the 
equivalence we observed exposes a hidden, unbroken supersymmetry of oo-bosons with 
super symmetric partners living on lattices of different sizes. 

5 Discussion 

This article is the third in a series studying the quantum mechanics of a simple super- 
symmetric matrix model. Designed originally to illustrate the usefulness of the large- 
approximation directly in terms of a Hamiltonian and a Hilbert space, the model turned 
out to have a very rich physics by itself, as amply illustrated in 

Here we have uncovered an even more intriguing aspect of this model: its connection, 
in the infinite 't Hooft-coupling limit, to one-dimensional statistical systems in which 
supersymmetry, if present, is very well concealed. We have found that our supersymmetric 
planar model is exactly equivalent to two such systems: the quantum XXZ Heisenberg 
chain, and a lattice gas of g-bosons. 

The intriguing pattern of strong coupling vacua discovered in the matrix model finds 
its explanation in the unusually simple ground states of the XXZ chain found by Baxter 
more than thirty years ago. Vice versa, the XXZ chain turns out to have a hidden su- 
persymmetry, which explains a host of degeneracies between seemingly unrelated energy 
eigenstates. Interestingly, the supersymmetry transformations change the number of lat- 
tice sites, so that different members of the supermultiplets live on different lattices. But 
even within a single sector/lattice, supersymmetry may turn out to be a powerful tool 
for studying the properties of the ground state and for understanding the meaning of the 
RS conjectures j8]. The fact that such ground states are annihilated by two supercharge 
operators should imply distinct properties for them, while finding an operator with the 
right algebra would allow to generate the full "staircase" of ground states starting from 
the lowest and simplest ones. 
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The second system, a gas of g-bosons in the hmit of an infinite deformation parameter 
g, is equivalent to our matrix model to an even greater degree than the Heisenberg chain. 
While the equivalence with the XXZ model holds only for a subset of all sectors of the 
Hilbert space, the cxo-bosonic gas is equivalent in all sectors of conserved boson number 
and for all lattice sizes. To the best of our knowledge, the specific Hamiltonian of the 
bosonic gas was unsolved till now; but in view of our chain of equivalences, the system 
should turn out to be actually soluble (e.g. via the Bethe ansatz for the XXZ chain). 

Finally, exact solubility of the XXZ model directly implies that same property for 
our large- iV, super symmetric matrix model at infinite 't Hooft coupling and, therefore, in 
"three quarters" of its fermionic sectors (since we have found no correspondence for the 
even-i?, even-F sectors). There are also indications that this solubility can be extended 
to the whole strong-coupling phase using the technique developed in |4]. The reverse is 
also an interesting question: are there statistical systems that can be mapped into our 
matrix model at finite 't Hooft coupling? 

Recently relations between the field theories and spin chains have become the subject 
of much interest and excitement in connection with the AdS/CFT correspondence, which 
opens a new way of studying gauge theories [TU]- |24j^ . In particular, the mapping between 
the A/" = 4 supersymmetric Yang-Mills theory and the XXZ chain with A = 1/2 has been 
discussed by Belitsky et al. [20j|. It is perhaps not an accident that we find a similar 
gauge-spin relation in our much simpler model. 
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Note added in proofs 

After this paper had been submitted we were informed by Jan de Gier that the existence 
of a hidden supersymmetry in the XXZ spin chain had already been pointed out in Ref. 
[25] and further studied in later work (e.g. in [26]). It looks however that, while in those 
papers supersymmetry is realized non-linearly just in terms of fermionic variables, in our 
case it takes a simpler form in terms of an equal number of bosonic and fermionic degrees 
of freedom. We are grateful to Dr. de Gier for this information. 

^We recall that, so far, the easiest applications of the AdS/CFT correspondence are also made in the 
A ^ oo (i.e. supergravity) limit. 
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